A Macroscopic Scale Model of Bacterial Flagellar Bundling 
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Escherichia coli and other bacteria use rotating helical filaments to swim. Each cell typically 
has about four filaments, which bundle or disperse depending on the sense of motor rotation. To 
study the bundling process, we built a macroscopic scale model consisting of stepper-motor-driven 
polymer helices in a tank filled with a high- viscosity silicone oil. The Reynolds number, the ratio of 
viscous to elastic stresses, and the helix geometry of our experimental model approximately match 
the corresponding quantities of the full scale E. coli cells. We analyze digital video images of the 
rotating helices to show that the initial rate of bundling is proportional to the motor frequency 
and is independent of the characteristic relaxation time of the filament. We also determine which 
combinations of helix handedness and sense of motor rotation lead to bundling. 

PACS numbers: 87.19.St, 87.16. Qp, 87.16.-b, 47.15. Gf 
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Introduction. Escherichia coli and Salmonella ty- 
phimurium cells have several helical propellers, or flag- 
ella, which they use to swim. Each flagellum consists of a 
rotary motor embedded in the cell wall, a short (50 nm) 
flexible hook which acts as a universal joint, and a heli- 
cal filament about 20 nm in diameter and about 10 /im 
long The trajectory of an individual swimming cell 
consists of runs interrupted by tumbles. For most of a 
run, the motors turn counterclockwise when viewed from 
outside the cell, the filaments wrap into a tight bundle, 
and the cell swims along a roughly straight path. Near 
the end of a run, one or more of the motors reverses direc- 
tion, the corresponding filaments come out of the bundle, 
and the cell moves erratically, or tumbles. The tumbling 
process is complex and involves polymorphic transitions 
of the filament first from the left-handed "normal" state 
to the right-handed "semi-coiled" state, and then to the 
"curly- P' state 0- The first transition reorients the cell 
body. When the motors resume their counterclockwise 
rotation, the curly- 1 filaments transform directly to the 
normal state and rejoin the bundle, and the cell resumes 
its initial speed 

The chemotaxis strategy of E. coli is to decrease the 
likelihood of tumbling during runs which happen to carry 
the cell toward higher concentrations of chemoattrac- 
tants. Thus, the formation and dispersal of the helical 
bundle is central to bacterial chemotaxis. Since the ra- 
dius of the flagellar filament is well below optical wave- 
lengths, and the motor rotation is relatively rapid (100 
Hz) , it is difficult to study the mechanics of the bundling 
process directly. Therefore, we built a macroscopic scale- 
model system consisting of flexible rotating helices in a 
very viscous fluid. By including the viscous fluid and 
properly accounting for the relative strengths of viscous 
and elastic stresses, our scale-model builds upon and ex- 
tends the work of Macnab, who studied the geometry of 
rotating flexible helices in a bundle 

Our paper begins with a discussion of the material pa- 
rameters of bacterial flagella and how we chose the pa- 
rameters for the experimental model. The next section 
describes the geometry of bundled helices and the sym- 



metry requirements for bundling, including helix handed- 
ness, motor-rotation sense and relative speed, and phase 
relations. Finally, we describe measurements of the char- 
acteristic time scale governing the initial stages of the 
bundling process. 

Scale Model. E. coli usually has several filaments per 
cell, but for the purpose of studying bundling it is sim- 
plest to consider the case with two filaments. Unless 
the two helices are identical and placed on diametrically 
opposite sides of the body (an unlikely occurrence), the 
hydrodynamic torque on the rotating helices causes the 
cell body to counter-rotate to make the total torque on 
the cell vanish. This counter-rotation plays some role 
in bundle formation 0, since it tends to wrap the 
filaments around each other. However, this mechanism 
for bundling does not lead to tight bundles @. Further- 
more, bundles readily form from filaments of stationary 
cells bonded to a surface Therefore, we disregard the 
effects of body rotation and translation and consider two 
helices rotated by stationary motors. We also neglect the 
effect of the cell body itself on the flow. 

At the characteristic length and time scales of bac- 
teria, viscous effects dominate inertial effects 0- To 
see why, recall that the Reynolds number Re = pv£/rj 
determines the relative importance of inertia to viscous 
stresses, where p is the density of the fluid, rj is the vis- 
cosity, v is a characteristic flow velocity, and £ is a charac- 
teristic length over which the flow varies. Using 30 /im/s 
as a typical swimming speed, 1 fim as a typical cell-body 
size, and the density and viscosity of water (p w 1000 
kg/m 3 , n « 10~ 3 N-s/m 2 ) yields a Reynolds number of 
Re w 3 x 10 -5 . (Using the typical helix radius I w 0.2 
/zm and the linear azimuthal helix velocity v « 2-7T x 100 
Hzx£ yields approximately the same Reynolds number.) 
At low Reynolds number, the flow field induced by a 
point force falls off inversely with distance @ , leading to 
a long-range hydrodynamic interaction between two ro- 
tating helices. (Since the net force on a swimming bac- 
terium vanishes, the far field flow induced by a swimming 
bacterium falls off faster than 1/r, roughly like a dipole, 
1/r 2 .) Another characteristic of low Reynolds number 



FIG. 1: Bundling sequence. The top panels show our left-handed model "semi- coiled" helices (see Table Pi at t = 0, 96, 168, 
and 264 s; the bottom panels show the helices at the same times, viewed from the side. The scale bars are 100 mm long; the 
helices are 310 mm long (from chuck to tip), 4.0 mm in diameter, and turning at 0.1 Hz. 



Form 


Pitch, /im 


Diameter, /im 


handedness 


normal 


2.3 


0.4 


left 


semi-coiled 


1.1 


0.5 


right 


curly 1 


1.0 


0.3 


right 


model normal 


11 x 10 4 


2.54 x 10 4 


left 


model "semi-coiled" 


6.4 x 10 4 


2.54 x 10 4 


right and left 



TABLE I: Helix parameters: approximate pitch, approximate 
diameter, and handedness for some polymorphic forms of flag- 
ellar filaments from and We also include the dimen- 
sions of our model helices. For the models, the diameter is 
the diameter of the mandrel. 



flow is that the drag per unit length on a long slender 
body depends weakly on the filament diameter and the 
shape of the cross section @- 

Elastic stresses balance the viscous stresses on rotat- 
ing bacterial flagella, causing the filaments to bend and 
twist as they wrap into a bundle. Filaments can also 
undergo polymorphic transformations when the motors 
turn clockwise; for simplicity, we do not attempt to in- 
clude polymorphic transformations in our macroscopic 
model. The bending modulus A = EI determines the 
resistance of an elastic rod to bending, where E is the 
Young's modulus and / is the moment of inertia of the 
cross section There are few measurements of the 

bending resistance, with reported values ranging from 



1(T 24 N-m 2 [13 to 1CT 22 N-m 2 pj. Likewise, the value 
of twist modulus C of the filament is not precisely charac- 
terized, although the twist compliance of approximately 
0.2 /im-long filaments attached to the hook and a locked 
motor has been measured |l3j . Since it is difficult to sep- 
arate the twist modulus of the hook from the modulus of 
the filament in this measurement, we will simply assume 
that C ~ A for the filament. For a filament with axial 
length L, the product of the characteristic filament relax- 
ation time and the motor angular velocity, M = r/ajL^/A, 
determines the importance of viscous drag relative to 
flexibility (see Q and jT^). Since A = EI oc Ea 4 , 
where a is the filament radius, M depends sensitively on 
the aspect ratio L/a. To estimate M for the bacteria, 
we will use L = 7 /im, the typical length of the filaments 
of reference 2J. Using uj = 100 x 27r rad/s, r\ = 10~ 3 
N-s/m 2 , and A = 10~ 23 N-m 2 yields M w 150 for the 
bacterial filaments. Note that the persistence length 0] 
£p = A/(k-Q,T) ks 2.5 mm is large, allowing us to disre- 
gard thermal fluctuations. 

There are three dimensionless numbers which deter- 
mine the helix shape: P/L, R/L, and a/L, where P 
is the helical pitch, and R is the helical radius. Note 
that the ratio R/P determines the pitch angle a via 
tana = lixRj P, and that P and L determine the number 
of turns (which need not be an integer). If we consider 
two helices with initially parallel axes, then the final ge- 
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ometrical parameter of interest is the spacing h between 
the two motors. 

Methods. Figure ^ shows the macroscopic model. Two 
stationary stepper motors under independent control ro- 
tate flexible polymer helices in a tank of silicone oil. The 
viscosity of the silicone oil is 100 N-s/m 2 , or about 10 5 
times the viscosity of water. The motors typically rotate 
at about 0.1 Hz, and the helix radius is about 13 mm. 
For these parameters, the Reynolds number Re ~ 10~ 3 , 
not quite as low as 10~ 5 but low enough to justify our 
neglect of inertial effects. To minimize the effects of the 
walls of the tank on the flow near the helices, all experi- 
ments are carried out with the helices near the center of 
the tank. The base of the tank is 420 mm by 420 mm, and 
the depth of the silicone oil is 330 mm. The helices are 
made by wrapping a Tygon™ tube around a cylindrical 
mandrel and filling the tube with epoxy. Food coloring 
is added to provide contrast in the video images. Once 
the epoxy hardens, the tube maintains its helical shape 
but is also flexible. We vary the flexibility by varying the 
radius of the Tygon™ tube. To get appropriate values 
of M in our model, we use helices with a filament radius 
of a few millimeters. The bending stiffness of each helix 
is approximately determined by measuring the deflection 
under a known force of a straight tube filled with hard- 
ened epoxy; for 2a — 4.0 mm, 4.7 mm, and 5.6 mm, we 
found A = (3.46±0.14) x 10~ 3 N-m 2 , (10.53±0.14) x 10~ 3 
N-m 2 , and (22.50 ±0.14) x 10 -3 N-m 2 , respectively. As 
in the case of the bacterial filaments, we assume C ~ A. 
Since the aspect ratio L /a for our model is much smaller 
than that of the full-scale bacterial flagclla, our model is 
not a true scale model. However, since the drag of a slen- 
der body depends only weakly on the filament radius , 
it is not necessary to attain exact geometric similarity. 
We used two representative helical geometries for our fil- 
aments (see Table [I] for a comparison of the dimensions 
of our model helices and bacterial flagella). We system- 
atically varied the stiffness of the "semi-coiled" model. 
Note that the (inner) radii of our model helices are all the 
same, since both types of helices were formed from the 
same mandrel, whereas the helical radius of the bacterial 
filament varies from polymorph to polymorph. All exper- 
iments reported here were carried out with a motor shaft- 
to-shaft separation of h = 61 mm, which would roughly 
correspond to the typical spacing of one micron in the full 
scale bacteria. To allay concerns about effects of fatigue 
and plastic deformation of the helices during the experi- 
ments, we compared the deflection of the helices under a 
known weight before and after some representative trials 
and found that these effects were insignificant. 

Bundling in the Scale Model. Figure shows a se- 
quence of snapshots of two rotating left-handed helices 
(with four turns each) at M 170 (see Movie 1, which is 
published as supporting information on the PNAS web- 
site, www.pnas.org). The induced flows cause large de- 
flections, and a bundle forms when the motors turn coun- 
terclockwise. Note that the helices wrap around each 
other in a right-handed sense; the flow field generated by 




FIG. 2: (a): Sequence of left-handed helices turning clockwise 
at 0.075 Hz, at times t = 102 s, t = 119 s, and t = 146 s. The 
helices cross, and the crossing point migrates up to the motor, 
where the helices periodically distort to escape a jam. 



each helix tilts the other helix, causing the helices to roll 
around each other and form a right-handed wrapping. 
In contrast with the case of two helices, the deflection 
of a single rotating helix with M « 170 is small, with 
the axial compression or extension less than ten percent 
of the axial length. But an understanding of the flow 
induced by a single helix sheds light on the flows that 
cause bundling. Consider a left-handed helix. As the 
helix rotates counterclockwise as viewed from the helix- 
side of the motor, apparent helical waves travel along the 
filament away from the motor and push fluid elements 
in the same direction. These fluid elements also rotate 
about the vertical axis of the helix in the same sense 
as the helix, counterclockwise, leading to a right-handed 
trajectory for each fluid element. 

Early scale-model experiments by Macnab using steel 
helices in air showed that left-handed helices twisted 
around each other in a right-handed manner can rotate 
indefinitely without jamming, but that left-handed wrap- 
ping leads to jamming for clockwise rotation and un- 
winding for counterclockwise rotation [3J. Macnab also 
pointed out that the state of lowest elastic stress of two 
left-handed helices with right-handed intertwisting is one 
in which the helices are coaxial, in phase, and in contact. 
Since this interwound state is unstable in the absence of 
external torque, Macnab used a guide at the end of the 
bundle to keep the helices from unravelling. In our exper- 
iment, viscous drag provides the stabilizing torques that 
keep the helices interwound. Note that the helices ro- 
tate against each other, leading to a region of high shear 
between the filaments. When the motors turn counter- 
clockwise, we find that the bundle persists indefinitely 
as long as the motor speeds are sufficiently low. If the 
motor speed is greater than 0.1 Hz, then the Y-shaped 
junction in the top right panel of Fig. ^ migrates up to 
the motors, causing a jam. Jams can also form when the 
motors turning the helices in a steady-state bundle re- 
verse; if a jam fails to form, the filaments simply unwind. 

In the example shown in Fig. ^ the two helices have 
no initial phase difference. Since the motors run at a 
common constant velocity and do not slip, the phase dif- 
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FIG. 3: Trajectories of (projected) helix tips in the a;y-plane. 
The spacing of the dots on each curve corresponds to one 
revolution. 

ference between the motor shafts is constant. If there 
is an initial phase difference, then the helices must bend 
and twist to attain the in-phase interwound state. Typ- 
ically one helix bends more than the other, leading to a 
bundle which is tilted relative to the axis of rotation of 
each motor. The motor characteristic in the scale model, 
constant speed, differs from that of the bacterial motor, 
which for typical loads and temperatures is roughly con- 
stant torque 0,0]. Thus, the phase- locking observed 
in the bacterial flagella may arise from a different mech- 
anism than that observed here. 

In the scale model, it is crucial for the rotation speeds 
of the two motors to be close for bundling to occur. If 
the motors speeds differ by as little as 8%, left-handed 
helices rotated counter-clockwise do not bundle. The he- 
lices cross each other, with the faster helix bending more 
than the slower helix. The axis of the faster helix waves 
back and forth, but never wraps around the slower helix. 

Since left-handed helices turning counter-clockwise 
bundle, symmetry implies that right-handed helices turn- 
ing clockwise will also bundle. Indeed, bundles of semi- 
coiled or curly flagella have been observed in fluorescently 
labeled filaments |2| . On the other hand, left-handed he- 
lices turning clockwise do not bundle. When the motors 
driving a bundle reverse, the bundle unwinds. Once the 
bundle has unwound, or if the (left-handed) helices are 
initially stationary, then the flow induced by the clock- 
wise rotation causes the helices to roll around each other 
and wrap in a left-handed sense. However, the wrapping 
is not nearly as tight as the counter-clockwise sense. The 
helices cross, and the crossing point migrates up to the 
motors; see Fig. |5J When the crossing point reaches the 
motors, instead of jamming, the helices distort at the 
proximal ends to allow the rotation to continue. None 
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TABLE II: Rotation speeds and ratios for the curves (labeled 
I- VII) shown in Fig. fU The columns labeled W\, u>2, and 
u>3 correspond to the 4.0 mm, 4.7 mm, and 5.6 mm tubes 
respectively. The empty spaces arise because most of the 
overlapping sets of curves (except for case IV) consist of just 
two curves. 




FIG. 4: Side view of (projected) tip-separation as a function 
of time for rotation speeds and bending moduli; see Table [2] 
for values. On each curve, there is a black dot for every 
revolution. The gold curves are for the 4.0 mm tube, the 
blue are for the 4.7 mm tube and the green are for the 5.6 
mm tube. 



of the other combinations of handedness and sense of 
rotation led to bundling in our scale-model experiment. 
Helices of the same handedness but turning in opposite 
directions tilt away from each other and also tilt away 
from the plane containing the initial helix axes. Helices 
of opposite handedness turning in the same sense cross 
and do not bundle; helices of opposite handedness turn- 
ing in opposite senses tilt away from each other. 
Characteristic Time Scale. There are two natural 
time scales governing the dynamics of our system: the 
motor period 27t/oj, and the elastic relaxation time scale 
77L 4 1 A. To determine which combination of these con- 
trols the initial rate of bundling, we measured the deflec- 
tion of each helix as a function of time. More precisely, 
we measured the position of the projection of the tip of 
each helix onto the local axis of rotation for the corre- 
sponding helix. Figure shows the trajectories of these 
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points projected to the icy-plane. In contrast to the case 
of Fig. we used helices with three turns [L — 240 mm) 
since with short helices it is easier to attain small val- 
ues of M. For M < 70, the viscous stresses induced by 
the rotating helices are too weak to deflect the helices 
much. The helices cross when M rj 70, and begin to 
wrap around each other when M rj 130. Since we fo- 
cus on the early stages of bundling, we do not show the 
complete trajectories for moderate to large M (M > 50). 
Defining A as the distance between ends of the helices in 
the side view (j/z-plane, where z is parallel to the initial 
helix axes, and x is parallel to the line connecting the 
two motors; Fig.^), dimensional analysis implies 

- = f(tot, M, Re, a/L, R/L, P/L, h/L). (1) 

Since we approximate C = A, we omit A/C from the 
list of dimensionless groups. In this paper we do not re- 
port how / depends on the geometrical parameters R/L, 
P/L, and h/L; we present results for semi-coiled helices 
(P/R RJ 5) with three turns, h = 61 mm, and L — 240 
mm. We varied rotation speed to and filament radius 
a. The function / depends on to implicitly through M 
and Re, and explicitly through tot. Since Re <C 1, we 
can approximate Re ri and consider / to depend on 
to only through M and tot. Likewise, / depends on a 
implicitly through M and explicitly through a/L. The 
explicit et/L-dependence of / arises from viscous drag be- 
cause the Young's modulus E = ira 4 /A enters through M. 
Since the drag on a slender body depends weakly on a/L, 
and since M depends sensitively on a/L, we consider / 
to depend on a only through M. Thus, for the fixed ge- 
ometrical parameters described above, we measure / as 
a function of tot and M. 

Figure 0] shows A as a function of tot for various rota- 
tion speeds and bending moduli, and Table ITU shows the 
speed corresponding to each A vs. tot curve. The ratio 
of the speeds for the overlapping gold and blue curves is 
roughly constant (see Table ITT|l . Note that the gold and 
blue curves labeled V do not overlap as closely as those 
labeled I, II, and III; correspondingly, lo\/uoi for the pair 
V noticeably differs from the ratios for I, II, and III. Like- 
wise, the blue and green curves labeled IV do not overlap 
as closely as those labeled III, leading to the discrepancy 
between the u>2 / to$ values for pairs III and IV. Although 
M systematically increases as the curve label varies from 
I to VII, there is a significant spread in the values of M 
for overlapping curves. This spread implies that our iden- 
tification of the bending modulus of a straight tube with 
that of a helix with the same diameter tube leads to an 



uncertainty greater than that of our measurements of A. 
Such an uncertainty is not unexpected since we neglect 
the effects of the helix twist modulus C for simplicity. De- 
spite the uncertainty in the helix bending modulus, the 
data of Fig. 0] and Table UJ] are strong evidence for the 
scaling behavior of the bundling helices: two helices of 
different stiffness but the same shape will yield the same 
A vs. tot curves for the proper ratio of rotation speeds. 
In fact, if we assume that the closely overlapping curves 
of Fig. 01 have the same M, then we can get a better esti- 
mate of elasticity ratios using A1/A2 = to±/t02, since we 
can measure the rotation speeds to high accuracy. 

A characteristic time of the early stages of the bundling 
process is the time at which the separation A reaches 
a maximum. Figure 01 reveals that the position of this 
maximum is proportional to the motor period, and de- 
pends only weakly on the elastic relaxation time scale. 
Thus, in the early stages of bundling, the helices follow 
the induced flow and the motor period determines the 
characteristic time scale. The helices begin to bundle af- 
ter about six revolutions, in rough agreement with the 
corresponding time for bacterial flagella Q • 
Conclusion. Our macroscopic scale model demon- 
strates that the bundling of bacterial flagella is a purely 
mechanical phenomenon, arising from the interplay of hy- 
drodynamic interactions, bending and twisting elasticity, 
and geometry. Our model allows us to study how the 
bundling phenomenon is affected by parameters which 
are difficult to control in the full-scale bacteria, such as 
the rate and direction of motor rotation. For fixed values 
of the geometrical parameters, we quantitatively char- 
acterized the early stages of the bundling process, and 
found that the initial rate of bundling is mainly deter- 
mined by the motor period. Future work should address 
the role of twisting vs. bending stiffness, the effect on 
bundling of the translational and rotational flow fields 
induced by a swimming bacterium, and bundle forma- 
tion when there are more than two filaments. 
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